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except the group of university students (P = 
0.028), were in Hardy-Weinberg equilibrium 
(by a M arkov chain method in GENEPOP). 
Furthermore, we compared the genotypes 
among individuals with SARS with differ¬ 
ent prognoses. If L-SIGN homozygosity is 
a protective factor against infection, it may 
also be associated with better prognosis after 
acquiring the infection. Therefore, we also 
examined whether homozygotes had a bet- 
m ter prognosis by classifying individuals with 
'■= SARS who had an uneventful recovery versus 
<E thosewho had severe disease and wereadmit- 
o ted to the intensive care unit for mechanical 
■3 ventilation support (an approach similar 
5 to that reported previously 5 ). However, we 

1 did not detect any significant association 
<d (P = 0.9, Supplementary Table 1 online). 

3 Sample size is the main limitation of both 
q studies. However, these two samples already 
| represent the few large' collections of indi- 
^ vidualswith SARS availablefor genetic study. 
S- To estimate the size of an overall effect, we 
^ performed a m eta-analysis of the two datasets 
§- together by the M antel-Haenszel test using 

2 controlgroupsinHardy-Weinbergequilibrium 
15 (twogroupsofcontrolsinthisstudyand ran- 
c dom controlsfrom Chan et al.; total n = 1,468; 

462 affected individuals and 1,006 controls). 

]5 The combined odds ratio was not significant 

if (combined OR = 0.84; 95% confidenceinter- 

£ val: 0.66-1.06, P =0.14). 

5 The difference in the results between the 

2 two studies was basically accounted for by 

o 
o 

CM 


a difference in the homozygote proportions 
in the controls (45.4% in this study versus 
55.0% in Chan et al.), while the homozy¬ 
gote proportions among individuals with 
SARS are almost identical (46.9% here ver- 
sus46.3% in Chan et al.).Thereason for the 
discrepancy in the homozygote proportions 
in the'control' groups is not clear. However, 
a subpopulation difference in allelic and 
genotypic frequencies exists between north¬ 
ern and southern Chinese. The seven-repeat 
allele was more prevalent in the Beijing 
sample (0.7 in Beijing versus 0.64 in Hong 
Kong; P =0.05), which also largely accounted 
for the higher proportion of homozygotes 
(55.7% in Beijing versus 46.0% in Hong 
Kong; P = 0.02). Unrecognized subpopula¬ 
tion structure may confound genetic asso¬ 
ciation studies. Results in thestudy by Chan 
etal. suggested that this confounding factor 
might be present. There were three groups 
of controls, includingtwo groupsof hospital 
controls (health care workers who worked in 
SARS wardsand affected individualsattend- 
ing various outpatient clinics) and a group of 
blood donor controls. I nteresti ngly, genotype 
distributions from both groupsof hospital- 
based controls deviated significantly (P < 
0.0001) or marginally (P =0.05) from Hardy- 
Weinberg equilibrium. 

In addition, other yet-unknown mecha¬ 
nisms (such as alternative splicing of the neck 
region, which could interfere with formation 
of homotetramersamong homozygotes) may 


account for the discrepancy between thetwo 
studies. Replication isan importantapproach 
to verify any significant genetic association 
findings 6 ' 7 , and additional association studies 
are required to establish the putative protec¬ 
tive effect of L-SIGN homozygosity against 
SARS or other infections. 
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Lack of support for an association between CLEC4M 
^homozygosity and protection against SARS 
coronavirus infection 


To the Editor: 

I n thej anuary 2006 issue, C han et al. 1 reported 
a significant association between severe acute 
respiratory syndrome (SARS) and a variable 
numberoftandem repeats(VNTR) polymor¬ 
phism in exon 4 of CLEC4M in a collection 
of individuals from Hong Kong. CLEC4M 
encodes L-SIGN ('liver/lymph nodespecific 
ICAM -3grabbing nonintegrin'), which serves 
as a receptor for many viruses, including SARS 
coronavirus (CoV) 2 .1 ndividuals homozygous 
forCLEC4M tandem repeats were reported to 
belesssusceptibletoSARSCoV infection.The 
authors also showed that cells homozygousfor 
CLEC4M repeats had a higher binding capac¬ 
ity for SARS CoV, higher proteasomedepen- 
dentviral degradation and alower capacity for 


trans infection. Thus, both genetic and func¬ 
tional studies suggested that homozygosity 
for CLEC4M was associated with protection 
against SARS CoV infection. 

Itisimportantto bear in mind that associa¬ 
tion studiesrequirereplication in independent 
populations 3 . We thereto re attempted to repli¬ 
cate thefindings of Chan etal. bygenotyping 
theVNTR polymorphism in three additional 
collections of case-control samples from 
northern China: (i) the'Beijing community 
population', consisting of 339 individuals with 
SARS and 227 random controlsrecruited from 
the community 4 ; (ii) the'Beijing health care 
worker (HCW) population', consisting of 42 
health care workers infected with SARSduring 
the course of hospital duty and 40 health care 


workers who had worked in SARS wards but 
remained freeof disease and wereconfirmed to 
be seronegative for SARS 5 and (iii) the'Tianjin 
population', consisting of 60 individuals with 
SARS and 129 disease-free controls (including 
85 random controls and 44 health care work¬ 
ers) 6 . Thethreecollectionsof case-control sam¬ 
ples an d th ei r ascertai n m en t cr i teri a h ave been 
described in detail previously (Supplementary 
Methods online) 4-6 . All groups except the 
individuals with SARS from the Beijing com¬ 
munity were in Hardy-Weinbergequilibrium. 
We found no significant differences in allele, 
genotype and homozygote or heterozygote 
frequencies between affected individuals and 
controls in thethree populations (Table land 
Supplementary Tablel online). Early reports 
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Table 1 Homozygote or heterozygote frequencies of the CLEC4M polymorphism in sam¬ 
ples from northern China 


VNTR polymorphism 3 


H eterozygotes 


H omozygotes 


Beijing community population 
All cases (n = 339) 

Random controls (n = 227) 

OR (95% c.i.) 

P value 


139 (41.0%) 
110 (48.5%) 


200 (59.0%) 
117 (51.5%) 


1.29 (0.89-1.87) 
0.19 


Cases with comorbid conditions (n = 52) 
Random controls (n = 227) 

OR (95% c.i.) 

P value 


20 (38.5%) 
110 (48.5%) 


32 (61.5%) 
117 (51.5%) 


1.36 (0.60-3.07) 
0.46 


Cases without comorbid conditions (n = 287) 
Random controls (n = 227) 

OR (95% c.i.) 

P value 


119 (41.5%) 
110 (48.5%) 


168 (58.5%) 
117 (51.5%) 


1.26 (0.86-1.85) 
0.23 


Severe cases (n = 19) 
Mild cases (n = 268) 
OR (95% c.i.) 

P value 


8 (42.1%) 
111 (41.4%) 


11 (57.9%) 
157 (58.6%) 


1. 00 (0.38-2.59) 
1 . 00 


Beijing HCW population 
HCWcases (n =42) 
HCW controls (n =40) 
OR (95% C.i.) 

P value 


14 (33.3%) 
18 (45.0%) 


28 (66.7%) 
22 (55.0%) 


1.51 (0.58-3.99) 
0.40 




Tianjin population 
All cases (n = 60) 

All controls (n = 129) 
OR (95% c.i.) 

P value 


33 (55.0%) 
72 (55.8%) 


27 (45.0%) 
57 (44.2%) 


1.18 (0.60-2.34) 
0.63 


|OR, odds ratio; c.i., confidence interval; VNTR, variable number tandem repeat; HCW, health care worker. 
a The heterozygotes are used as the reference group, and all ORs and P values are adjusted for age and gender. 
Primer sequences used for genotyping are listed in Supplementary Table 2 online. 


have shown that some nongen etic factors, such 
as comorbid conditions (including diabetes 
mellitus, hypertension, heart disease, tuber¬ 
culosis, asthma and malignancy), are risk fac¬ 
tors for the development of SARS 7 ' 8 and may 
confound the contribution of genetic factors 
to thisdisorder. H owever, after stratification by 
comorbid conditions, theassociation remained 
nonsignificant in our Beijing community 
population (Table 1). Of the287 affected indi¬ 
viduals without comorbid conditions, 19 were 
individuals with severeSARS who wereidenti- 
fied by their admission to intensive care units 
or by their death, and the remaining 268 were 
individuals with mildly symptomatic SARS. 
To account for this, we assessed whether there 
was an association between homozygosity for 


CLEC4M tandem repeatsand severity of SARS, 
but we found none (Table 1). This result may 
be due to the limited number of individuals 
with severeSARS in thecurrent study and will 
require confirmation in additional studies. 

Collectively, theresultsin ourthreecollections 
of case-control samples from northern China 
arenotsupportiveofthefindingsofsignificant 
association between theVNTR polymorphism 
and SARS risk reported by Chan et al. 1 . There 
are several possible reasons for the inconsis¬ 
tent results. First, inadequate power may bean 
explanation of our negative results. However, the 
samplesizein our Beijing community popula¬ 
tion (with 339 affected individualsand 227 con¬ 
trols) is approximately similar to that used by 
Chan et al. (with 285 affected individualsand 


380 random controls), and thissamplesizehad 
power >0.94 to repl icate the effects by C han et 
al. (calculated bythegenetic power calculator 9 ). 
Additionally, just by using our Beijing commu¬ 
nity sampleset, wehavesuccessfully confirmed 
thepositiveassociations between themannose- 
binding lectin polymorphisms and SARS risk 4 
that wereobserved previously in a case-control 
sample from Hong Kong 10 . Furthermore, the 
consistency of the negative associations in our 
Beijing H CW population and our Tianjin popu¬ 
lation strengthensour results. 

Second, there may be a small, population- 
specific difference in the contribution of 
CLEC4M polymorphism to SARS susceptibil¬ 
ity. This might occur if there were population 
differences in linkage disequilibrium pattern 
or allele frequencies of CLEC4M . Indeed, the 
homozygote or heterozygote frequencies in 
individuals with SARS in our Beijing com¬ 
munity or Beijing HCW populations were 
significantly different from those reported for 
affected individuals in Hong Kong community 
or Hong Kong HCW populations (P =0.0066 
and P = 0.017, respectively; Supplementary 
Table 1). Additionally, the difference in allele 
frequency between Hong Kong outpatient 
controls and Tianjin random controls and 
the differences in homozygote or heterozy¬ 
gote frequencies between Hong Kong HCW 
controls and Tianjin HCW controls were also 
significant (P =0.0064 and P =0.017, respec¬ 
tively; SupplementaryTablel). Furthermore, 
there was also a significant difference in the 
genotype and homozygote or heterozygote 
frequency between the Chinese population 
and Europeans 1 . 

Another possibility isthatChan etal. unwit¬ 
tingly neglected to account for the potential 
confounding factors that may distort thecon- 
tribution of CLEC4M VNTR polymorphism to 
SARS susceptibility. Although Chan etal. took 
stringent precautionsto stratify their samples 
by 'health care worker', population samples 
may still differ by many other factors that 
depend on setting and context of recruitment, 
such as age of the subjects at SARS onset, sex 
of the subjects and any comorbid conditions. 
Unfortunately,Chan etal. do notprovidedata 
with regard to such important information on 
confoundingfactors. Last, theinitial findingsof 
Chan etal. may not represent real associations 
and might befalse positives. In genetic associa¬ 
tion studiesof common diseases, thereisa very 
low prior probability of detecting a true asso¬ 
ciation result when accounting for statistical 
adjustmentfor multiplecomparisons. I ndeed, 
the genetic associations presented by Chan 
et al. were marginally statistically significant 
(P = 0.027, 0.045 and 0.031, when compar¬ 
ing all SARS samples to random controls, 
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community SARS to outpatient controls and 
HCW SARS to HCW controls, respectively). 

Thus, wedid not find any significant differ¬ 
ences in allele, genotype and homozygote or 
heterozygote frequencies between cases and 
controlsin our three independent populations 
of northern Chinese. Although the biological 
plausibility of L-SIGN and thefunctional evi- 
denceof theVNTR polymorphism in theorigi- 
nal report remain interesting, weurge that the 
association between CLEC4M polymorphism 
and SARS be investigated in other subpopu¬ 
lations of ethnic Chinese origin (for example, 
TaiwaneseorGuangdongChinese) or in those 
of different ancestry, such as Europeans. 
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Chan etal. reply: 

Tang etal. and Zhietal. report that in 
contrast to our findings 1 , they were unable 
to find association between homozygosity 
or heterozygosity of the CL EC 4M (L-SIGN) 
exon 4 tandem repeat polymorphism and 
SARSCoV susceptibility in the Chinese 
population. Their data cannot conclusively 
negate our findings for the reasons below. 

We agree with Tang et al. that the 
difference between their analysisand ours 
is largely accou nted for by the difference 
in the percentage of homozygosity in 
controls. Tang’s controls included neonatal 
cord blood, healthy elderly individuals 
jZT»v aged >70 years, local university students 
pJand control samples collected in Beijing. 

T he homozygosity and heterozygosity 
distribution of their Beijing controls is very 
similar to that of our Hong Kong random 
controls. H owever, the former three control 
groups collected in Hong Kong are poorly 
matched by age for comparison with 
individuals with SARS and our random 
controls (Supplem entary Table 1 online), 
and such an age distribution clearly excludes 
the population aged 25- 69 years, which 
makes up the largest proportion (>70%) 
of SARS-infected individuals during the 
HongKongoutbreakin 2003 (ref. 2). Age 
and gender are well-known confounding 
factors in any case-control study, and for 
SARS, individualsaged >65 years and <18 
years are actually associated with a lower 
risk 3 .Thus, it issurprising that Tang et 
al. chose these specific control groups at 
extremeage ranges, given that there are no 
obviou s obstacles for recru iti ng appropriate 


age-matched controls in HongKong.They 
also failed to perform logistic regression 
accounting for age differences, which should 
have been incorporated in their statistical 
analysis. T he age and gender distri bution 
of the SARS and control populations in our 
study, in contrast, is more properly matched, 
and by logistic regression, our results 
remain statistically significant (Table 1 and 
Supplementary Table 1). 0 ur results remain 
consistent after accou nting for comorbid 
conditions (Table 1). 

Tang etal. also used control groups 
differingwidely in age to determineif age 
had any effect on genotype frequencies. 

H owever, thei r data per se al ready showed 
a significant difference in overall genotype 
distribution, when comparing neonates 
versuselderly individuals versus university 
students(P = 0.029, /} test). T he difference 
in genotype between neonates and university 
students is also significant (P = 0.009 
by CLU M P), suggestingthat age- related 
selection may exist for CLEC4M genotypes. 
Indeed, age-dependent variation of allele 
and genotype frequencies has been reported 
for other genes 4 ’ 5 . 

Tang etal. also claim to show a 
subpopulation difference in allele and 
genotype frequencies between northern 
and southern C hi nese. T hey note that "the 
seven-repeat allele was more prevalent in 
the Beijing sample (0.7 in Beijing versus 
0.64 in H ong Kong, P = 0.05), which also 
largely accounted for the higher proportion 
of homozygotes (55.7% in Beijing versus 
46.0% in H ong Kong, P = 0.02)." It seems 
these quoted Hong Kong figures refer to 
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frequencies from their cord blood group 
alone. Such a comparison isdifficultto 
justify,again because age isseriously 
mismatched. It should be noted thatthereis 
no difference in either the seven-repeat allele 
frequency or the proportion of homozygotes 
between the H ong Kong random controls 
of Chan etal. 1 and theBeijingcontrolsof 
Tang etal. (SupplementaryTable2 online). 
Barreiro etal. 6 ’ 7 investigated the CLEC4M 
homozygote and heterozygote distribution 
of different ethnic groups and reported 
that the proportion of homozygotesin 
East Asians, consisting predominantly of 
Chinese, was53% (SupplementaryTable 
3 online), afiguresimilar to the Hong 
Kong random controls of C han et al. and 
theBeijingcontrolsof Tang etal. There is 
also no significant difference in allele or 
genotype frequencies of these two control 
populations from theEast Asian samples 
of Barreiro etal. 7 (SupplementaryTable 
2). W hen these resu Its are taken together, 
little subpopulation structure is observed 
forCLEC4M homozygote and heterozygote 
distribution between the northern and 
southern Chinese populations. 

F i nai ly, given the apparent absence of 
subpopulation structure for CLEC4M 
homozygote and heterozygote distribution 
in the C hinese population, we performed 
a m eta-analysis of our data set and that 
of Tang et al. by the M antel-H aenszel test 
usingall control groupsthatarein Hardy- 
Weinbergequilibrium: our random controls 
and outpatient controls (n = 670), and 
cord blood, healthy elderly individuals 
and Beijing controls(n = 827) of Tang 
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